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Abstract 

Maximal 't Hooft loops are studied in SO (3) lattice gauge theory at finite temperature T. 
Tunneling barriers among twist sectors causing loss of ergodicity for local update algorithms are 
overcome through parallel tempering, enabling us to measure the vortex free energy F and to 
identify a deconfinement transition at some f3^ lt . The behavior of F below f3^ lt shows however 
striking differences with what is expected from discretizations in the fundamental representation. 
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Topology plays an important role in the non-perturbative dynamics of Yang-Mills theo- 
ries. In particular the vacuum condensation of topological excitations might explain quark 
confinement and the existence of a mass gap through various scenarios described in the 
literature [H 0, HI . Some of these models allow to define a topological order parameter for 
the finite T deconfinement transition; for 't Hooft magnetic vortices, classified by the first 
homotopy class Z N of the continuum Yang-Mills gauge group SU(N)/Z N the change in 
free energy F = AU — TAS for their creation might play such role and has received broad 
attention, in particular in lattice discretizations at zero and finite T 0, 11 E[. The main 
problem in such non-perturbative regularizations is that creating a vortex is equivalent to 
the introduction of a non-trivial twist. For discretizations in the fundamental representation, 
transforming under the enlarged gauge group SU(N) (N = 2 in this paper), this cannot be 
implemented dynamically but only via a modification of their boundary conditions (b.c), 
the generalized partition function Z being defined through the weighted sum of partition 
functions with fixed twisted b.c. Since each of them must be determined by independent 
simulations their relative weights can only be calculated through indirect means 

Universality arguments are often cited to claim that results in lattice Yang-Mills theories 
will not depend on the discretization chosen. A natural alternative in calculating F would 
therefore be to directly discretize the theory in the SU(N)/Zn representation (SO (3) in 
our case) with periodic b.c, naturally transforming under the continuum Yang-Mills gauge 
group (7[. The adjoint partition function Z((3a) should in fact be equivalent to Z provided 
the 5*0(3) native constraint a c = F|p g9c sign(Tr F t/p) = 1 is satisfied for every elementary 
3-cube c, where Up denotes the plaquettes belonging to its surface dc 0, [§]: Z 2 magnetic 
monopoles are suppressed and only closed Z 2 magnetic vortices winding around the bound- 
aries are allowed, i.e. in this limit adjoint actions dynamically allow all topological sectors 
which in the fundamental case must be fixed through twisted b.c. Moreover, since the 
standard spontaneous center symmetry breaking argument for the deconfinement transition 
does not apply to the center-blind adjoint discretization the question whether F behaves 
as an order parameter is of major interest, also in light of recent studies for alternative 
descriptions of confinement in centerless theories 0. 

A practical obstacle one needs to overcome in investigating the adjoint theory is the 
appearance of a bulk transition at some /3a, separating a strong coupling chaotic phase (I) 
continuously connected with the fundamental action, where (a c ) ~ 0, from a weak coupling 
ordered phase (II) where (a c ) ~ 1 (loL 11, I n phase II, where one wishes to exploit 



the relation between Z((3a) and Z mentioned above, high potential barriers separating twist 
sectors suppress tunneling among them for local update algorithms [7]. On the other hand 
a well-known result is that a center blind Z 2 monopole suppression term in the action 
A 5^ c (l — er c ) weakens the order of the bulk transition while moving it down into the strong 
coupling region (see the curved dashed line in Fig. QJi) [12]. For asymmetric lattice sizes 
N T x Ng, N T <C N s indications for a finite temperature critical line f3^ lt (X, N T ) within phase 
II (horizontal dashed line in Fig. for N T = 4) have already been found from simulations 
at fixed twist H H " 



In this paper we solve the problem of ergodic updates through parallel tempering (PT) 
[l7 . 18l | for the adjoint Wilson action with Z 2 monopole suppression 

S = (3 A J2 (l - ^ Tr ^) + A X> " *0 > C 1 ) 
where Up denotes the standard plaquette variable and Tr^O = (Tr F 0) 2 — 1 = Tr^(0 2 ) + 1 
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FIG. 1: (a) Phase structure in the (3 A - A plane with paths of couplings indicated as used in our 
simulations, (b) Adjoint Polyakov loop average (La) versus 0a at fixed A = 0.8 and N T = 4. 



the adjoint trace. Our PT paths in the (3a — A plane extend over the second order leg of 
the bulk transition straight into the inner region of phase II as shown in Fig. with two 
paths at fixed Pa, one path with fixed A, and one path consisting of a fixed-/?^ and a fixed- 
A part. In this way we avoid the necessity to cross high potential barriers. This enables 
us to account for all twist sectors, restoring ergodicity, and to study F at finite T. As 
a qualitative indicator of deconfinement we have also measured the adjoint Polyakov loop 
(La) = C^2 s Tr A L(x))/(3Ng), L(x) = Ylt U&(x, t), which strictly speaking cannot behave 
as an order parameter for the deconfinement transition, delivering at most information 
on the screening length for the effective adjoint potential (see Fig. |IJ)). A preliminary 
report on the present investigation was presented in 3| ] further results on other observables 



as well as a detailed description of the implementation of the algorithm will appear in a 
separate publication j^. To our knowledge this is the first successful attempt to study 
the deconfinement transition of the 5*0(3) lattice gauge theory via an ergodic simulation at 
large volume as well as the first determination of F in the confined phase of a center-blind 
discretization of Yang-Mills theories. 

In PT |T7|,ll8( we update K configurations i — 1, . . . , K, with couplings (A, p A )i swap- 
ping neighboring pairs from smaller to larger values according to a Metropolis acceptance 
probability, satisfying detailed balance, P S wa,p(hj) = min[l, exp(— ASy)], where 

ASij = S[(X, p A )i, Fi] + S[(X, p A )j,Fj) ~ S[(X, p A )i, T 3 \ - S[(X, p A )j, H ■ ( 2 ) 



Compared to other methods |21|, |22| the appeal of PT is its easy implementation both at 
criticality and away from it, needing no knowledge of re- weighting factors or other dynamical 
input. This is a welcome property for us since we wish to go as far as possible from the 
bulk transition into phase II. In view of various experiences with simulated tempering one 
also expects PT to be more efficient than multi-canonical simulations. For the success of 
the method in the case under consideration the softening of the bulk transition to 2 nd order 
is crucial, since we "transport" tunneling from phase I at lower A into phase II at larger 
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A, where twist sectors are well defined but frozen. To work at low A, i.e. through a 1 st 
order bulk transition, would make barriers too high and kill any hope of ergodicity at large 
volume, as experienced in 0] for iV s > 8. 

Some care is of course necessary also with our method. In particular to maintain a suffi- 
cient swapping acceptance rate u the distance between neighboring couplings must diminish 
with the volume. On the other hand to keep cross-correlations under control one does not 
wish the acceptance rate to be too high. We have chosen to tune the parameters for each 
path and volume at hand so to keep the acceptance rate roughly fixed at around to = 12%, a 
value for which we empirically find a good balance between auto- (in the sense of freezing of 
the sectors) and cross-correlations. For the paths in Fig^i contributing to the left and right 
branches in Fig [2K details on the statistics, i.e. number of configurations N and number of 
ensembles K, are given in Table |U The paths at fixed (3a = 0.4 and 0.65 for N s = 16 shown 
in Fig |2b were calculated with K = 7 and iV = 30000 configurations. To remain on the 
safe side, in Fig|2K we have chosen not to quote the ensembles related to the end points of 
the paths since, having no further configurations to swap with, they might be affected by 
systematic errors. We wish however to stress that this has never been observed in the PT 
literature and we also have no indication that this might be the case. Further details on 
the algorithm and a detailed analysis of correlations will be reported in a forthcoming paper 
|2o| . For the fixed A paths of Fig. [TJi, along which main simulations have been performed, 
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TABLE I: Statistics achieved for the data given in Fig |2k. The left branches correspond to (a) 
while the right branches together with pieces at fixed (3a = 0.95 and varying A correspond to (b). 

we can fit very well the step 5(3a needed to keep u fixed with a law of the form 

5(3a(^N s )^^-, (3) 

s 

where a(12%) = 2.15(3) in the (3a = 0.95 — 1.09 range considered, although we expect it to 
change with the (A, (3a) window. Such scaling implies that in order to explore a fixed region 
A(3a of parameter space the number of ensembles will scale like 

K ~ ^N 2 S . (4) 
a 

Therefore we cannot go too deep into phase II since the number of PT simulations will 
eventually become too large for the computational means at our disposal. 

The native 50(3) temporal twists are given by Zi = N~ 2 J2 Xj , Xk Yip e plane i,4 s ig n Trp£/p, 
(eijfc4 = 1) for i — x,y,z @, Q|- Within phase II they are well defined having values close 
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to either 1 or -1 for each configuration. The partition functions_restricted to fixed twist are 
given as expectation values (i>k) = Z\ z= k/Z of the projectors 



i ~ n [i + sign(^)] , v x = i Yi n t 1 + ( x - 2< m si s n (^)] > 

i=x,y,z j=x,y,z i=x,y,z 

l E II [1 - (1 - 2«S«) fflgn(zi) , ^3 = ^ II [1 - sign(^)] ■ (5) 



J=x,y,z i=x,y,Z 



t=x,y,z 



The 2- and 3- twist sectors can of course only exist on T 3 |3j, ly] . From Eq. , since for an 
adjoint theory a change of twist sector leaves the action unchanged, AU = and: 



F — —T log 



3Z n 



log 



3(^o> ' 



(6) 



The factor three in the denominator is again due to the three equivalent 1-twist sectors on 
T 3 rather than one as on R 3 0,0]. With such a choice F will be zero if all twists are equally 
probable, i.e. on top of the bulk transition and everywhere in phase I, where twist sectors 
are however ill-defined due to the presence of open vortices. Eq. (jUJ) obviously implies F = 
in the T — > limit 0, 0, S] as long as the AS contribution remains bounded. Fig. |2t shows 
numerical results for F/T in lattice units at fixed A = 0.8 obtained along the two separated 
paths indicated in Fig. QJi. Errors are given combining statistical errors with auto and cross 
correlations. The data start on top of the 2 nd order bulk transition and go upward to what 
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FIG. 2: (a) aN T F versus Pa at A = 0.8 and various N s . (b) aN T F versus A — A c at fixed (5a for 
N s = 16. X c (Pa) denotes the position of the bulk transition. 



we interpret as the finite T deconfinement transition - given the rapid growth of ( L 



15, 16] . 



Fig. Eh) as well as of F and taking fixed- twist results for varying N T into account 

The behavior in phase I and just at the bulk transition as well as in the deconfined 
phase at large (3a (upper phase II) is in agreement with the standard vortex arguments 
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for confinement if vortices behave "chaotically" then F should be zero and the theory 
confines, while as deconfinement occurs F should rise as F ~ aN%, cr(T) being the dual 
string tension. This is qualitatively in agreement with our data. While F ~ close to the 
bulk transition, we find a strong rise of F for (3a > (3 A nt = 1.01(1) allowing to locate the 
finite-temperature transition in phase II. As already explained above, we are not able to 
go too deep into the deconfined phase with the computing facilities at our disposal. Thus, 
we cannot really check whether the data for (3 a 3> (3 A nt are consistent with the expected 
O(Ng) plateaus or whether they get saturated with respect to the thermodynamic limit, 
i.e. whether a can indeed be calculated. The effort necessary to this purpose, assuming 
that the estimate in Eq. (j3J) still works at higher (3a and even taking into account that for 
higher volumes the asymptotic behavior should kick in earlier, we would need to simulate 
around 50 parallel ensembles for each volume, again for a statistics of at least O(10 5 ) per 
configuration in each ensemble. For volumes with N s > 20, for which finite size effects start 
to be reasonably small, this goes beyond the computational power at our disposal, although 
it should be manageable with a medium sized PC cluster. 

For (3a < (3 C A lt i i- e - throughout the confining region of phase II, we find negative values 
for F which stabilize at large 3- volumes (N s > 20). F < comes as a surprise, meaning that 
vortex production is enhanced as compared to phase I. This is in contrast to what expected 
from arguments valid within the fundamental representation, i.e. F(T) = throughout the 
confined phase. For an independent check we have carried out simulations at a few (fixed) 
lower /^-values, varying A i.e. at somewhat lower temperatures (see the horizontal paths 
drawn in Fig. |T^l). The results are plotted in Fig. |2b- Again we find F < 0, but the plateau 
values do not behave monotonously as a function of (3a- Passing some minimum they increase 
again for decreasing (3a- This is compatible with the expectation that the free energy should 
go to zero in the zero-temperature limit 0, 0, i| . A systematic extrapolation for different 
volumes and N T would be required to confirm this behavior and to decide whether F/T 
itself vanishes or goes to a constant value. 

Let us draw the conclusions. The main result of the present paper is the success in 
sampling the full partition function via ergodic PT Monte Carlo simulations and determining 
the free energy F for the creation of a Z 2 vortex in pure SO (3) Yang-Mills theory at finite 
T. We have seen a clear indication for a deconfinement transition consistent with earlier 
findings of a second order transition at fixed twist. Furthermore we find that F does not 
vanish in the confined phase at T ^ 0, vortex creation being enhanced throughout it; F 
cannot therefore serve as an order parameter in a strict sense. This implies that the adjoint 
theory is unable to exhibit an order parameter for center symmetry breaking in any form, 
much like in the case of strictly centerless groups |9j. This is not in contradiction with the 
confining properties of the model, F = being a sufficient but not necessary condition for 
confinement away from the T = limit 0, • Moreover while vortex suppression for T < T c 
would have been difficult to justify in light of the literature 2J] the vortex enhancement we 
observe does not contradict that they can play a role in describing confinement, although 
one cannot speak of vortex condensation in the usual understanding. 

Let us finish with a short remark on the universality problem we seem to face. When 
identifying the partition functions Z ~ Z((3a), invoking universality for observables, one 
should be cautious. First of all, although Z 2 monopoles become suppressed in the con- 
tinuum limit of Z, (cr<-) is still far from unity for the range of parameters commonly used 
in the simulations jj, y, E| and open vortices might still dominate the partition function. 
Moreover Z((3 A ) simply does not allow to define physical observables in the fundamental 
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representation. Expectation values of fundamental Wilson and Polyakov loops and all their 
correlators vanish identically for all /3a, i-e. a fundamental string tension cannot be defined 
in a straightforward way. To our knowledge bounds for F from reflection positivity or the 
connection with the electric flux have only been derived within the fundamental represen- 
tation 0,0, An interpretation of universality implying that any observable will assume 
the same value in any discretization is therefore trivially contradicted by the above consid- 
erations. However, a slightly more conservative reading can agree with our findings: the 
truly physical properties measurable in "experiments" like glueball masses and the critical 
exponents at the transition should be reflected by physical observables which can be defined 
irrespective of the discretization chosen. 

We thank Ph. de Forcrand, J. Greensite, E.-M. Ilgenfritz, T. Kovacs, M. Pepe and U. J. 
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